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' Abstract 

This paper generalizes both the binary Deutsch-Jozsa and Grover algorithms to n- 
^ , valued logic using the quantum Fourier transform. Our extended Deutsch-Jozsa algorithm 

is not only able to distinguish between constant and balanced Boolean functions in a single 
1/^ ' query, but can also find closed expressions for classes of afRne logical functions in quantum 

oracles, accurate to a constant term. Furthermore, our multi- valued extension of the Grover 
algorithm for quantum database search requires fewer qudits and hence a substantially 
smaller memory register, as well as fewer wasted information states, to implement. We 
note several applications of these algorithms and their advantages over the binary cases. 



§ ! 1 Introduction 

Quantum computers have the ability to implement algorithms in inherently multi-valued logic. 
Multi- valued quantum logic gates have already been developed by Muthukrishnan and Stroud 
, [T], and shown to be feasible in a linear ion trap model. Generalizing quantum algorithms to 

^ ■ multi-valued logic can also present marked advantages over other binary classical and quantum 

! algorithms. 

m 

. The original binary Deutsch-Jozsa algorithm 2 considers a Boolean function of the form 

/ : {0, 1}'' {0, 1} implemented in a black box circuit, or oracle, Uf. Input states are put in 
a quantum superposition as query (x) and answer (y) registers so that their state vectors are 
' expressed in terms of the dual basis 0] 

^' 

|0') = ^(|0) + |1)) and |1') = ^(|0) - |1)), also denoted |+) and |-) 
V 2 V 2 



o 



X 



\^ • The oracle is defined by its action on the registers: Uf\xy) — \x)\y © /(a;)), where the |a;) 

5^ \ register is the tensor product of input states |a;i) • • • \xr)- When it is promised that the function 

in question is either constant (returning a fixed value) or balanced (returning outputs equally 
among and 1), the algorithm decides deterministically which type it is with a single oracle 
query as opposed to the 2''"^ -j- 1 required classically. The corresponding circuit is shown in 
Figure 1 {/"^ denotes r wires in parallel). 

Grover's algorithm [5] is a probabilistic search algorithm usually presented in the context of 
searching an unsorted database. It consists of the iteration of a compound "Grover operator," 
which consists of an oracle and an inversion-about-the-average operator D, on a superposed 
register of search states as well as an ancillary qubit. The oracle marks the searched state by 
rotating its phase by tt, while the inversion-about-the-average operator increases the amplitude 
of this state. 
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Figure 1: The Deutsch-Jozsa circuit 

In this paper, we present extensions of both the Deutsch-Jozsa and Grover algorithms to 
arbitrary radices of multi-valued quantum logic. We denote addition over the additive group 
Z„ by the symbol and the Kronecker tensor product by (g). 
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Figure 2: The Grover circuit 

The Hadamard transform is a special case of the quantum Fourier transform (QFT) in 
Hilbert space T-Ln- The well-known Chrestenson gate for ternary quantum computing is also 
equivalent to the Fourier transform in "Ha- Cereceda [3] has generalized the Deutsch algo- 
rithm using two qudits to d-dimensional quantum systems where d = 2^. However, placing no 
restrictions on the number of computational basis states n leads to a far more versatile charac- 
terization of the Deutsch-Jozsa algorithm. The Fourier matrix of order n over the n^^ roots of 
unity uj'^ = e*^^'''/" is as follows: 



F = 

7), 



Definition 1 The action of the quantum Fourier transform is described by QFT^ : \j) ^ 
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\k) for j G [4]. The QFT matrix can be expressed as 



^ n— 1 n—1 



\j){k\. 



2 The Deutsch-Jozsa Algorithm in Radix n 

We will use the QFT in place of the Hadamard transform in the multi- valued equivalent of the 
Deutsch-Jozsa algorithm. First, we define the functions that we are working with. 
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Definition 2 An r-qudit multi-valued function of the form 

/:{0,l,...,n-ir^{0,l,...,n-l} 

is constant when f{x) = f{y) Vx, y € {0, 1, . . . , n — 1}'' and is balanced when an equal number 
of the domain values, namely n''^^, is mapped to each of the n elements in the codomain. 

In multi- valued logic, there are n constant functions mapping each element in Z„ to a fixed 
element and n! balanced permutative (bijective) mappings of single-qudit inputs. For functions 
on r qudits, there are accordingly n''!/n'"~^!" balanced (surjective) mappings. 

Theorem 1 All affine functions defined as f(xi, . . . ,Xr) = Aq (B Ai with 
Aq,. . . ,Ar e Z„ are either constant or balanced functions of r qudits. 

Proof. Those affine fimctions for which all coefficients A^^o = are constant. For affine 
functions with at least one nonzero coefficient of Xi, the domain is {0, 1, . . . , n — 1}'", which can 
be identified with {0, 1, . . . ,n'" — 1} = IJ^ 5m where Sm = {m,n+m,2n+m, . . . ,n'^ — n+m,} for 
every m G {0, 1, . . . , n — I}. Each Sm is a set of size n'^'~^, each of whose elements reduces modulo 
n to m. Since f{p) = f{q) if p = q (mod n), every element in the codomain {0, 1, . . . , n — 1} 
is assigned to exactly n^~^ different elements in the domain. Such affine functions satisfy the 
definition of a balanced function. 



The proof of the n-ary Deutsch- Jozsa algorithm will be aided by a simple lemma: 

n-l 

Lemma 1 Primitive n*^ roots of unity w satisfy w"'^ = for integers a ^ (mod n) . 

fe=o 

Proof. Consider the polynomial z" — 1 = 0, of which is a root. 1 being a real root for all 
integers n, this can be factorized as {z — l){z"~-^ + z""^ + • • • + 1) = 0. Since a = r (mod n) 
for some r where < r < n — 1, = co^ ^ 1. Therefore, X^^Zg z'^ = for z = co" . 



This leads to our main result. 



Theorem 2 The n-ary Deutsch- Jozsa algorithm applied to multi-valued functions of r qudits 
can both distinguish between constant and balanced functions with a single oracle query and 
determine a closed expression for an affine function in Uf, excepting the constant term, as 
follows: 

1. The constant term Aq is preserved in the phase of the x-register at output {liJ~^°), which 
is lost during measurement. 

2. The coefficients Ai,...,Ar are determined by the state of the x-register at output, 

\Ai,...,Ar). 




n-l) 



tlV'l) tlV'2) t|^3> tlV'4) 

Figure 3: The n-ary Deutsch-Jozsa circuit 
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In practice, the y-register would not be measured, but we follow through with the relevant 

calculations to demonstrate that its state at the output is constant, regardless of the function 
in the oracle. The x- and y-registers are written separately as factors of the entire tensored 
state of the circuit at each step 

First, we consider the case in which the function f{x) hidden in the oracle is constant: 

ivi) = io)«ni) 

® +1 1 1 n—1 

^ n^ — 1 n — 1 

^ 1^3) = ^ E 1^) ® ^ E ® /(^))- 

" x=0 * y=0 

At this point, we can transfer the action oiUf from the basis states themselves onto their phases 
by observing that if a basis vector \j ® k) is appended with the phase factor ^ , then \j) itself 
must have phase cp^^'^. This yields: 

-J — 1 ^ n—1 

^ n^— 1 n—1 

x=0 ^ y=0 

Because wc assume our function to be constant, e-*2ir/(a;)/Ti ^^^^ regarded as a global phase 
factor. Subsequently, the QFT on the x-register can be computed explicitly: 

n^—l n^ — 1 

Fn^-^ = ^J2T. giving 

v"' j=0 fc=0 

IV'3) IV'4) = 

j=0 ;c=0 a;=0 j=0 k=0 y=0 

In the standard basis, the inner product {k\z) = dkz- We can hence reduce the above to 

^ n'^ — ln^ — l ^ n— In— 1 

1^^) = J_e-i27r/(x)/„ ^ ^ ei2xife/„|^.^ ^ i^^ei2,rOei)fc/„|^.^_ 

j=0 k=0 ^ j=0 k=0 

By Lemma 1, all basis states \j) in the .x-rcgister will have null amplitudes for j ^ 0. Similarly, 
all basis states \j) in the y-register will have null amplitudes for j ^ n — 1. It follows that 

|^4) = |0)«nn-l) 
with a phase factor of e~*2ir/(a;)/ra £qj. constant functions f{x). 

The balanced case is similar. After initializing and superposing our states as above, we 
obtain: 

— 1 -| n—1 

IV3) = e-^2,r/(x)/n. X ^ ^ ^i27^y/n<y 



x=0 ^ y=0 
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In this case, the phase factor e~'^'^/(^V" cannot be assumed to be global because its value is 

dependent \ipon x. The outpiit of the r/-register will be the same as in the constant case, so we 
need only proceed with the state of the x-register. After applying the second QFT: 

® 1 n^—l 

j=0 k^O x=0 



j=0 a;=0 

It is now necessary to show that jx — f(x) = some constant C\ or f{x) = jx — C, for a fixed 
value of j 7^ and all x in the domain {0, 1, . . . , n'' — 1}. This would allow e*^'^'-^/", or uj'^' , 
to be the phase factor of some basis state \j) other than jO)'^'' as in the constant case, with a 
unit probability of measurement. Equivalently, since addition and multiplication are modular, 
/(xi, . . . , Xr) = —C ® jiXi © • • • ® jrXr must hold for some j G {1, n'' — 1}, ji representing 
the i*"^ digit of j. This is the definition of a non-constant affine function, so j exists only when 
the hidden function in the oracle is affine. The x-register is therefore measured to be: 

f{x) being balanced for j 7^ and constant otherwise. In consequence, the Deutsch-Jozsa 
algorithm gives a deterministic output only when f{x) is restricted to being either constant or 
balanced, and affine (we will give an example of the algorithm applied to a non-affine function 
below). However, observe that the state jO)'^'' will have a zero probability of measurement for 
all balanced functions, whether affine or not, because its amplitude is 

2;— a;— 2:— 

by Definition 2 and Lemma 1. Thus, the algorithm is still deterministic in the sense that it can 
always distinguish between constant and either affine or non-affine balanced functions, although 
with no fixed output in the latter case. 



Finally, the generalized Deutsch-Jozsa algorithm has a useful property that becomes appar- 
ent in multi- valued logic; it can not only distinguish between constant and balanced functions, 
but can determine explicitly the function j{x\, . . . , Xr) = ffi ^12:1 ® • • • ® A^Xr implemented 
by the oracle excepting the constant term ^0, given that it is affine. As calculated above, 
the constant term Aq of such a function is encoded in the phase of the x-register at output 
(^0 = —C, where the phase factor is oj*^), while the respective coefficients Ai,. . . ,Ar oi x are 
determined by the basis vector \j) = \Ai, . . . , Ar). Since the phase of the a;-register is lost at 
measurement, only Aq cannot be retrieved. We regard affine functions that differ only in the 
constant term as a "class." 



The y-register can be ignored for ease of theoretical computation (not in practice) if the 
oracle, corresponding to the diagonal operator 

directly encodes the action of f{x) into the phase of the a;-register. We use this scheme below. 



Example 1 (Deutsch-Jozsa for an affine function) 

Uf contains the following balanced function defined on two qutrits \AB): 
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We begin at IV'3), after the states have been initiahzed. 
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Lemma 1 is used for simphfication (1 + w + = 0): 
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=^ |V4)=C^'|21), 

from which we derive a closed expression for the afEne function /(xi, X2) = Aq® Aixi ® ^20:2 
in Uf, save for the constant term, by taking {2,1} as the respective coefficients Ai and A2 
(although theoretically, Aq should be (—2) mod 3 = 1): 

f{xi, X2) = (Constant) ® 2xi ® X2- 



Example 2 (Deutsch-Jozsa for a non-affine function) 

Uf contains the following balanced function defined on two qutrits \AB): 
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Again beginning at IV'3): 
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1^4)= 3|01) + 
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3 -102) + -111) + -112) + -121) + 3 

The basis state |12) hence has a 4/9 probabihty of measurement, with ah others having 1/9 
probabiUty. To determine that this balanced function is non-affine, enough measurements of 
the a;-register are required so as to obtain different states at output. Furthermore, observe that 
if the function is not affine, there is still a relatively high probability of measuring the output 
state of the algorithm that would be associated with a similar afhne function. For example, 
compare the Marquand chart of the non-afBne function in {// above with that of the afRne 
function f{xi,X2) = xi ® 2x2 associated with the output state |12): 
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Because an afhne multi- valued function in our context is defined in terms of the modulo- 
additive operator ©, any arbitrary function is afhne iff it satisfies the cyclic group property 
(this is easy to determine when the outputs are plotted in a Marquand chart as above, which 
requires that rows and columns be successive cyclic shifts of each other), and the state after Uj 
is factorable - e.g., entanglement does not occur in the Deutsch-Jozsa circuit. 

Our extended Deutsch-Jozsa algorithm requires only one measurement to deterministically 
distinguish an expression for an afHne function of radix n and r inputs up to the accuracy of a 
constant, given that it is affine. 



3 The Grover Algorithm in Radix n 

The quantum Fourier transform could also be applied to Grover's quantum search algorithm, 
which, in a multi- valued context, would reduce the number of qudits necessary in the memory 
register for large search spaces. Specifically, while a search space of N elements in the binary 
Grover algorithm would require [IgA''] qubits to implement, a radix-n circuit would require 
[log„ N] qudits. This reduces the number of qudits by a factor of about Ig n. 
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Take \'4)) = (J>i|0))'^'' as the vector representing the 7'-qudit search space of N elements 
in an equal superposition, N = n^, and |jo) as the searched state. The action of the oracle 
Uf\xy) = \x)\y © fix)) gives us reason to redefine the (now ill-termed) inversion-about-the- 
average, or diffusion, operator as 

D^{l-L,-')\^){i^\-I 
and the selective phase rotation operator U\ig) as 

I-{l-LU-')\t0){t0\, 

if only the action on the x-register is considered. It is easily verified that both of these operators 
are unitary. The expression (1 — ll>~^) appears because Uf, adding f{x) to the basis states y 
modulo n, accordingly rotates the phase of the searched state by — 27r/n. f/|i(,) can be con- 
structed from Uf hy transforming the j/-register, initialized to via QFT to the superposed 
state X]fc=o as shown below [TO] . 




fe=o ^ fc=0 

Figure 4: Implementation of U\ig^ 

Again, if the y-register is ignored for computational purposes, then the oracle takes the same 
form as U^i^^ : 

where f{x) — 1 if 2; — iq and f{x) = otherwise, thereby marking the searched state through 
phase rotation. 

Hence, the entire Grover iteration can be expressed as the unitary operator 

G = = (1 - io-')mw + h){io\) - ^i^^^i^)(zoi - /. 



TV 

For practical implementation, the multi-valued D operator can also be decomposed into a QFT, 
a selective phase rotation of the state |0) by — 27r/n, and an inverse QFT: 

(1 ~ - / = J-r[(l - c.-i)|0)(0| - = -Tf:U\,,^{Tl)^^ 

Since the multi-valued Grover algorithm deals with rotations through complex angles in 
n-dimensional Hilbert space whereas the binary algorithm deals only with rotations through 
real angles in a real subspace of 2''-dimensional Hilbert space, the matter of precisely how 
many times to iterate is more complicated. Although each Grover iteration does rotate \'^) 
through a complex angle Q where cos 6* = {^\G^^^ G^\'^) after fc iterations, the interpretation 
of such a rotation as increasing the component of \^) in the direction of |io) ([2], [10]) is not as 
straightforward in higher radices. It is easily seen, though, that as n becomes very large, the 
number of necessary iterations increases since the (1 — u)~^) factor in G approaches 0. Thus we 
find a tradeoff between time and space complexity when going to multi- valued logic. 

We now argue that our multi- valued extension of the Grover algorithm still converges to the 
desired state in 0{^/N) steps. In the subspace of 'H®^ spanned by (the non-orthogonal basis) 



8 



{|io), IV)}) we have the eigendecomposition 
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Suppose A'' ^ 1. We calculate the amplitude of |io) after k iterations. Keeping only terms of 
lowest order in 1/A'', namely 1/2, and assuming fc <C iV over the course of the evolution that 
we consider (this latter assumption will be justified at the end of the computation), we have 



(io|G*^|V) = ( 1 1/^ )M ( ^+ ) ( ? ) 



y-x 



1 + 



-uj — 



y-x 



-UJ 



a; -V^(l-a;-^) 
-V2(l-a;-i)^ 



,^-1/2(1 



TV 



^^N y/N 



l\k-l 



[k{l - LJ-^) - LJ-^] . 



All the terms in uj are of bounded and small magnitude; hence k comes to dominate fairly 
quickly, and we see that the amplitude of |io) in the state of the register grows linearly with 
k. Furthermore, since the coefficient of k is on the order of 1/VN, 0{V N) steps are required 
to increase the magnitude of the amplitude of the target state beyond any given threshold, as 
desired, consistent with the assumption that fc <C AT for the entire evolution. 



4 Applications and Conclusion 

We have demonstrated that the multi-valued Deutsch-Jozsa algorithm is able to test for the 
constancy or affinity of a multi-valued logical function in a single step, preserving exponential 
speedup over classical algorithms. Furthermore, the multi-valued Grover algorithm allows for 
a reduction in memory register size by using fewer qudits to represent more information states, 
reduces the number of "wasted" information states by searching unsorted databases of close to 
arbitrary size, and still affords a quadratic speedup over classical algorithms. 

While the original Deutsch-Jozsa algorithm is mainly of theoretical interest, its multi- valued 
extension could potentially find application in image processing to distinguish between maps of 
colored patterns represented by affine functions in a Marquand chart. A salient feature of texture 
recognition is identifying and comparing predictably repeating patterned images. Such images 
could be easily encoded by maps of affine functions with each pixel represented by an output 
value, and each color represented by a diff'erent radix-n digit, n corresponding to the number of 
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colors in the image. Since we have shown that the multi- valued Deutsch-Jozsa algorithm gives 
a unique deterministic output for any class of afRnc function (i.e., differing in only the constant 
term) in 0(1) time, it presents a very efficient method of characterizing a texture image by a 
unique "signature" (the output state) without using the substantial processing power required 
by classical means. Since textures repeat, the constant term would be insignificant since it serves 
only to permute the entries in the function map. The image would only have to be encoded 
in the oracle, which would require a quantum computer with programmable gate arrays, and 
could be realized with prospective technology. 

Also, Grover's algorithm is well-known for its pertinence to a wide range of NP problems for 
which it could be adapted by modifying the oracle architecture. Several practical problems that 
could be solved include graph coloring, the traveling salesman problem (to which the Grover 
algorithm could be applied to verify whether any path has a cost less than a given threshold), 
and the general NP-complete Boolean satisfiability (SAT) problem. For example, in searching 
for a valid coloring of a planar map (where no two adjacent regions are of the same color) from 
the set of all possible colorings, any coloring of the corresponding graph with s vertices could 
be encoded in a string of s radix-4 digits; hence, the 4-valued Grover algorithm would be most 
efficient in terms of memory register size. Grover's algorithm is very versatile in that it can be 
modified for any general class of search problem; given a memory register of all possible inputs, 
the oracle can check each element against a Boolean or multi- valued logical function constructed 
from a certain number of constraints. Hence, most applications of Grover's algorithm are based 
on the reducibility of many important problems to SAT. The SAT problem could potentially find 
even more application when reformulated for multi- valued logic to allow for different numbers of 
parameters, using modulo-n addition and multiplication as natural generalizations of Boolean 
XOR and AND. 

The fact that the number of iterations needed to effectively implement Grover's algorithm 

increases with radix size leads to the question of the "optimal" radix for Grover's algorithm. This 
is a topic for further research. However, from numerical analysis in MATLAB, we conjecture 
that ternary logic is optimal for quantum searching in that it requires the smallest number of 
iterations for a search space of any given size. This may be related to its closeness to base e, 
the most "natural" base for computing. 
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Appendix 

In the multi- valued Deutsch-Jozsa algorithm, for radices higher than two in which Uf is merely 
r + 1 wires in parallel (implementing the constant function / : {0,1}'' 0), the initialized 

state |0)®''|1) is mapped to a different output |0)'^''|n — 1). Practically, this makes no differ- 
ence because the output of the y-register is discarded. However, Theorem 3, and the related 
characterization of J^^, make clear why this is so. 

Theorem 3 Four iterations of the QFT gives an identity mapping (the binary Hadamard gate 
is a special case that is also self-inverse). 
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Proof. Let denote the q entry of Tn- Row p and column q of are given by 

n—1 n—1 

(p| = Y.u'^iP-^Hkl and |g) = ^u;'=(^-i)|fc), 

respectively. We have [^^]pg — {p\q) = ^ J2k=o 0;'^^^+'"^^; p + q — 2 is some integer a that is 
nonzero when (p+q) mod n ^ 2, so by Lemma 1, = 1 if (p + q) mod n = 2 and [^^]pg = 

otherwise. Then all such indices p,q G {1,2, ... ,n} for which [J^n]pq — 1 satisfy p + g = Con + 2 
for some Co € Z„. Co is further restricted to {0,1} because any larger values of Con + 2 exceed 
the maximum value of p + q, or 2n. Therefore, either p + q = 2oTp + q = n + 2, giving the 
solution sets 

p = (7 = 1 or {p, q} = {n — m + 1, m + 1} for m G {1, 2, . . . , n — 1}, 
which correspond to the permutation matrix 

n-1 

= \n — m)(m|, \n) taken to mean |0). 

m=0 

Thus 



^ j^O k^o m^O 
^ n — 1 n — 1 

^ n — 1 n — 1 

= ^EE^""""'"ij>(fci 



V j=0 k=0 

= n 

by our previous arguments, and given that Fn is unitary, 

F"^ = F = I 

is immediate. This is a property that the QFT shares with the continuous Fourier transform. 
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